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13.  AASTRACT 


The  dependence  of  the  second-  and  third-order  polarizabilities  on  ground-state  polarization 
was  measured  for  a  series  of.  donor-acceptor  polyenes  using  electric  field  induced  second 
harmonic  generation  and  third  harmonic  generation,  respectively.  The  changes  in  ground-state 
polarization,  associated  with  the  donor/acceptor  strength  or  solvent  polarity,  were  probed  by  x-ray 
crystallography,  nuclear  magnetic  resonance,  electronic  absorption  and  Raman  spectroscopies. 
The  observed  behavior  of  the  second-  and  third-order  polarizabilities  as  a  function  of  ground-state 
polarization  agrees  well  with  theoretical  predictions.  In  particular,  positive  and  negative  peaks,  as 
well  as  sign  changes,  were  observed  for  both  the  second-  and  third-order  polarizabilities.  The 
dependence's  for  the  second-  and  third-order  polarizabilities  are  consistent  with  a  derivative 
relationship  between  them.  In  addition,  the  third-order  polarizability  of  a  series  of  molecules 
possessing  zero  bond  length  alternation  was  found  to  be  negative,  in  agreement  with  predictions 
based  on  the  relationship  between  the  polarizabilities  and  ground-state  geometry. 
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1.  ANNUAL  TECHNICAL  REPORT 

a)  A  comprehensive  list  of  objectives  of  the  research  effort 

The  objectives  of  this  work  are  to  obtain  an  accurate  assessment  of  the  limitations  of  the  second 
hyperpolarizability  (y)  for  molecules  of  a  given  length,  to  develop  design  strategies  for  materials 
wherein  y  has  been  optimized  for  a  given  length  and  to  synthesize  molecules  with  optimized  y.  In 
order  to  do  so,  a  combined  theoretical  modeling,  synthesis  and  measurement  program  has  been 
initiated. 

b)  Status  of  the  research  effort. 

Recent  studies  have  suggested  that  optimizing  the  polarizabilities  of  donor-acceptor 
polyenes  requires  a  specific  degree  of  mixing  between  neutral  and  charge-separated  resonance 
structures.  The  degree  of  mixing  is  related  to  the  donor/acceptor  strength  and  can  be  correlated 
with  the  ground-state  polarization  (i.e.  the  dipole  moment)  and  a  molecular  parameter,  bond 
length  alternation  (BLA),  defined  as  the  difference  between  the  average  single  and  double  carbon- 
carbon  bond  lengths  in  the  polymethine  chain.  Figure  1  illustrates  how  chemical  modifications 
can  systematically  tune  the  bond  length  alternation  by  varying  the  degree  of  mixing  between 
neutral  and  charge-separated  structures.  For  imsubstituted  polyenes,  or  chromophores  with  weak 
donors/acceptors.  A,  the  ground-state  structure  resembles  that  of  the  neutral  resonance  form 
resulting  in  large  positive  BLA.3  Increasing  the  mixing  between  the  two  limiting  resonance  forms 
by  utilizing  the  stronger  dicyanovinyl  acceptor,  B,  results  in  increased  ground-state  polarization 
and  smaller  BLA.3  Complete  mixing  of  the  limiting  resonance  forms  results  in  essentially  zero 
BLA,  as  for  symmetrical  cyanines,  C.^  In  neutral  donor/acceptor  systems,  increasing  the  donor 
and/or  acceptor  strength  further,  D,  results  in  the  charge-separated  resonance  form  being  the 
dominant  contributor  to  the  ground  state,  resulting  in  negative  BLA.5 

For  a  particular  chroraophore,  the  degree  of  BLA  can  be  fine-tuned  by  utilizing  solvents  of 
varying  polarity  since  the  degree  of  mixing  between  neutral  and  charge-separated  resonance 
structures  is  known  to  be  sensitive  to  this  perturbation.3.6,7  since  polar  solvents  can  stabilize  the 
charge-separated  structure  to  a  larger  extent  than  nonpolar  solvents,  the  ground-state  polarization 
of  a  donor-acceptor  molecule  can  thus  be  increased  by  utilizing  increasingly  polar  solvents,  since 
this  increases  the  contribution  of  the  charge-separated  form  to  the  ground-state  geometry. 
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Figure  1.  Contribution  of  neutral  and  charge-separated  resonance  forms,  given  qualitatively  by 
the  area  of  the  open  and  shaded  circles  respectively,  to  the  ground  state.  BLA  values,  tuned  by 

varying  donor/acceptor  strengths,  were  determined  from  x-ray  diffraction.3-5  Et  h  C2H5. 

Quantum  mechanical  calculations  have  been  used  to  predict  the  dependence  of  the 
polarizability  of  donor-acceptor  molecules  on  bond  length  alternation.^^S  Xhe  predicted 
correlation  of  the  first-,  second-  and  third-order  polarizability  (a,  P  and  y  respectively)  with  bond 
length  alternation  is  depicted  in  Figure  2  where  it  can  be  seen  that:  (i)  the  polarizabilities  can  be 
optimized  as  a  function  of  BLA;  (ii)  relationships  between  the  various  order  polarizabilities  exist, 
which  are  consistent  with  the  n^h-order  polarizability  being  related  to  the  derivative  of  the  (n-l)th- 
order  polarizability.  In  this  paper,  the  polarizabilities  of  a  series  of  donor-acceptor  polyenes 
exhibiting  a  wide  range  of  ground-state  polarization  were  probed  to  test  the  structure-property 
relationships  depicted  in  Figure  2. 


THHORETICAT .  DEPENDENCE  OF  TFIE  POLARIZABILITIES  ON  BLA 


A  detailed  explanation  of  the  predicted  dependence  of  the  polarizabilities  on  BLA  is  presented 
elsewhere,^  but  briefly  it  can  be  understood  by  consideration  of  the  behavior  of  the  terms  in 
simplified  sum-over-states  expressions  for  a,  p  and  y,  as  a  function  of  BLA. 


(1) 


3 


where  |j.  and  E  are  the  dipole  matrix  element  and  transition  energy,  respectively,  between  the 
subscripted  states.  The  subscripts  g,  e  and  e'  label  the  ground,  first  excited  and  upper  excited 
states  respectively. 
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Figure  2.  Schematic  illustration  of  the  predicted  dependence  of  a,  P  and  y  on  bond  length 
alternation.  The  ground-state  polarization  increases  from  A  to  E. 
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The  first-order  polarizability,  a,  is  dependent  on  two  terms,  (l/Egg)  and  [igg,  which  both 

peak!  at  the  cyanine-limit  of  zero  bond  length  alternation  (Figure  3),  explaining  the  dependence  of 
a  on  BLA  in  Figure  2.  The  two-level  model^  for  P  is  similar  to  an  a-term  scaled  by  a  factor  of 
(^ee-Pgg)/Ege-  The  difference  in  dipole  moment  between  the  ground  and  excited  state,  (Hge- 
|igg),  exhibits  positive  and  negative  peaks  at  BLA  values  corresponding  to  regions  A  and  E 
respectively  (Figure  3).1»2,8,10  jhus  p,  as  a  function  of  increasing  ground-state  polarization, 
starts  positive,  increases,  peaks  in  a  positive  sense,  decreases,  crosses  through  zero  at  the  limit  of 
zero  BLA,  becomes  large  and  negative,  peaks  in  a  negative  sense  and  eventually  decreases  in 
magnitude  (Figure  3).1>2,8,10  p  displays  positive  and  negative  peaks  closer  to  the  zero  BLA  limit 
(i.e.  regions  B  and  D,  respectively)  than  the  (liee’i^gg)  (Figure  3)  since  (l/Egg)  and  Pge 

peak  at  the  zero  BLA  limit. 

“■(^A 


abode 


Bond  Length  Alternation 

Figure  3.  Dependence  of  (pige)^.  (l/Ege)^,  (|iee-Pgg)  and  p  on  the  ground-state  structure.  The 
ground-state  polarization  increases  from  A  to  E.  The  point  where  (iigg-iigg)  and  P  are  zero 

corresponds  to  the  bond-equivalent  cyanine-limit  of  zero  bond  length  alternation.  The  curves 
were  calculated^  for  Me2N-(CH=CH)4-CHO,  geometry-optimized  in  the  presence  of  a  static 

electric  field  (used  to  modify  the  geometry  from  the  highly  bond-alternate  polyene-limit  (region 
A),  through  the  bond-equivalent  cyanine-limit  (region  C)  to  the  bond-alternate  charge-separated 
polyene-limit  (region  E)),  at  the  semiempirical  INDO  level  using  a  sum-over-states  approach  in 
combination  with  an  INDO-configuration  interaction  technique  that  includes  full  single  and  limited 
double  configuration  excitations. 
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The  three-level  model^  for  the  third-order  polarizability,  y,  consists  of  a  negative,  a^-like 
term,  Nl,  a  two-photon  term,  PI,  and  a  P-like  term,  P2.  For  chromophores  with  large  positive 
BLA,  such  as  unsubstituted  polyenes,  the  P2  term  is  negligible  due  to  small  (liee'M-gg) 
term  dominates  Nl  resulting  in  positive  y  (Figure  By  increasing  the  ground-state 

polarization  and  thus  decreasing  BLA,  the  P2  term,  like  p  (Figures  2  and  3),  starts  to  increase  and 
hence  y  increases.  However,  increasing  the  ground-state  polarization  toward  the  cyanine-limit 
(region  C)  also  increases  INI  I  and  therefore  a  positive  peak  in  the  y  curve  is  evident  at  larger 
magnitudes  of  BLA  than  the  peak  in  the  p-like  P2  term.  Upon  further  polarization,  P2  peaks  in  a 
positive  sense  and  starts  to  decrease  and  IN II  continues  to  increase  and  thus  y  decreases  in 
magnitude  until  INII  =  PI  +  P2  and  y  is  zero.7>8  At  the  bond-equivalent  cyanine-limit  (region  C), 
the  P2  term  is  zero  since  (liee'ligg)  is  zero.  The  y  curve  exhibits  a  negative  peak  at  this  point 
since  both  INII  and  PI  peak  but  INII  >  IPll.7,8  To  a  first  approximation,  the  behavior  of  yin  the 
region  of  negative  BLA  mirrors  that  in  the  positive  BLA  region  (Figure  2). 


The  observation  that  a,  p  and  y  appear  to  be  related  by  a  derivative-like  relationship  can  be 
understood  by  considering  a  Taylor  series  expansion  of  the  dipole  moment. 


where 


li  =  lio  + 


0 


o 


(4) 

(5) 


so  that 


(6) 


Now  if  we  consider  a  molecule  to  be  polarized  by  a  strong  field,  F,  that  can  modify  the  degree  of 
bond  length  alternation  associated  with  the  molecule,  the  dipole  moment  will  depend  on  F  in 
some  characteristic  way, 

|X  =  ti(F)  (7) 

for  which  a  single  Taylor  series  with  constant  coefficients  is  not  a  sufficient  description. 
Qualitatively,  we  would  expect,  and  AM-1,^  INDO^  and  valence  bond  calculations^^  show:  (i) 
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that  for  weak  fields,  such  that  the  BLA  is  polyene-like,  n  will  be  small  and  weakly  dependent  on 
field;  (ii)  that  when  the  field  is  large,  such  that  the  structure  is  close  to  the  bond-alternate  charge- 
separated  polyene-limit,  p.  will  be  large  and  weakly  dependent  on  field;  (iii)  that  for  moderate 
fields,  |i  will  be  intermediate  and  will  vary  most  rapidly  with  field,  i.e.  the  rate  of  change  of  the 
BLA  with  field  is  largest  at  the  cyanine-limiL 


At  any  particular  value  of  F,  say  F',  we  can  expand  p.  in  a  Taylor  series,  and  thus 


a  = 


aF 


F' 


p  = 


aa 

aF 


ap 

aF 


F' 


(8) 


Therefore  the  functions  a(F),  p(F)  and  7(F)  will  be  given  by 


aF 


aF 


aF 


(9) 


leading  to  derivative  relationships  between  a  given  polarizability  and  the  lower-order 
polarizability.  Keeping  in  mind  that  the  ground-state  polarization,  and  hence  BLA,  is  a  function 
of  the  field  strength,  i.e.  BLA  =  BLA(F),  that  if  we  want  to  map  the  polarizabilities  as  a  function 
of  BLA,  that  the  change  of  variables  is  nonlinear,  i.e. 


9p  9F 

aF  aBLA 


(10) 


since  BLA  is  a  nonlinear  function  of  F.  This  leads  to  a  change  of  the  shape  of  the  a,  P  and  y 
curves  versus  BLA  compared  to  those  versus  F.  Nonetheless,  the  relation  of  the  peaks  and  zero 
crossings  of  the  different  polarizabilities,  e.g.  the  peak  of  a  and  the  zero  crossing  of  P,  are 
preserved  and  we  will  use  these  special  points  to  test  experimentally  the  validity  of  these 
relationships. 

MODEL  COMPOTTNDS  AND  DEGREES  OF  BLA 

The  molecules  (Figure  4)  studied  to  test  the  predicted  structure-property  relationships  were 
chosen  since  strong  evidence  exists  that  they  approximately  cover  the  BLA  range  A-E  in  Figure 
2.10  The  bond  length  alternation  values  were  coarsely  changed  by  utilizing  donors  and  acceptors 
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of  varying  strength  and  fine-tuned  using  solvent  polarity.  Bond  length  alternations,  determined 
from  x-ray  crystallography,  of  unsubstituted  polyenes  and  3,  containing  the  weak  formyl 


acceptor,  have  been  reported  previously^* and  allow  an  assignment  of  1-3  to  region  A.  This 
assignment  is  consistent  with  the  positive  solvatochromism  displayed  by  3  in  all  solvents  used 


(Table  I). 

1 


■< - ► 


Positive 

BLA 


2 


- ► 


Figure  4.  Canonical  resonance  structures  for  the  donor-acceptor  polyenes  investigated.  Electron 
donor/acceptor  strength  increases  from  3  to  8.  Et  =  C2H5  and  Bu  =  n-C4H9. 


Compound  4  with  the  stronger  dicyanovinyl  acceptor  possesses  a  bond  length  alternation  of 
-t-O.OlSA  in  the  solid  state,  assigning  this  molecule  to  region  B.  Interestingly,  this  small  BLA 
value  suggests  that  in  polar  solvents  the  ground-state  geometry  is  cyanine-like  whereas  in 
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nonpolar  solvents,  due  to  a  decreased  contribution  of  the  charge-separated  state,  the  ground-state 
structure  is  more  similar  to  an  unsubstituted  polyene.  These  solvent-dependent  geometry  changes 
on  a  related  molecule  (B,  Figure  1)  were  probed  by  non-resonant  Raman  spectroscopy  (Figure 
5).3  The  Raman  spectrum  of  the  unsubstituted  polyene,  (Figure  5a),  exhibits  two  main  bands, 
those  due  to  the  C-C  (1150  cm'l)  and  C=C  modes  (1600  crn'l).  The  Raman  spectrum  of  the 
bond-equivalent  cyanine,  (Figure  5d),  exhibits  increased  activity  at  intermediate  frequencies 
consistent  with  the  loss  of  discrete  C-C  and  C=C  character.  As  predicted,  the  Raman  spectrum  of 
the  dicyanovinyl  moiety  in  nonpolar  solvents  (CHClj^  Fijgnre  5b)  is  Qualitatively  similar  to  that  of 
the  polyene  whereas  the  spectrum  in  polar  solvents  (CH3CN,  Figure  5c)  is  qualitatively  similar 


1000  1200  1400  1600  1800  2000  1000  1200  1400  1600  1800  2000 
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Figure  5.  Raman  spectra  for  (a)  The  polyene  decatetraene  (1,  Figure  4)  in  the  solid  state;  (b)  B 
(Figure  1)  in  CHCI3;  (c)  B  (Figure  1)  in  CH3CN;  (d)  A  cyanine  (C,  Figure  1),  in  CH3CN.3 

The  counterion  for  the  cyanine  was  perchlorate. 

The  ground-state  polarization  can  be  increased  (bond  length  alternation  decreased)  further  by 
utilizing  the  barbituric  and  thiobarbituric  acid  acceptors,  5  and  6  respectively.  In  fact,  BLA  for  6 
in  the  solid  state  is  -0.014A  (Figure  1),  which  indicates  that  the  ground-state  geometry  is  biased 
slightly  towards  the  charge-separated  canonical  resonance  form.  Both  5  and  6  display  positive 
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solvatochromism  in  nonpolar  solvents,  whereas  negative  solvatochromism  is  displayed  for  5  in 
highly  polar  solvents  and  6  in  moderate  polarity  solvents,  consistent  with  5  being  in  region  B/C 
and  6  in  region  C. 

Molecules  7  and  8  display  negative  solvatochromism  in  all  solvents  used.  Furthermore,  the 
large  IH-^H  coupling  constants  (Figure  6)  across  the  central  carbon-carbon  bond  are  consistent 
with  a  trans  double  bond  (as  depicted  in  the  zwitterionic  form  in  Figure  4,  right)  assigning  7  and 
8  to  regions  C/D  and  D/E  respectively.  These  coupling  constants  increase  with  increasing  solvent 
polarity,  consistent  with  an  evolution  of  the  ground-state  geometry  towards  a  more  charge- 
separated  structure. 

The  product  of  the  ground-state  dipole  moment,  ji,  and  the  second-order  polarizability,  p, 
was  determined  by  solvent-dependent,  electric  field  induced  second  harmonic  generation  (EFISH) 
measurements  at  1907  nm  using  apparatus  and  methodology  described  elsewhere.  The  third- 
order  polarizability,  y,  was  determined  by  solvent-dependent  third  harmonic  generation  (THG) 
measurements  at  1907  nm.  A  description  of  the  experimental  apparatus  and  data  reduction 
scheme  has  been  reported  previously. 


#  barbituric  acid,  7. 
thiobarbituric  acid,  8. 


Solvent  E^(30) 


Figure  6.  The  IH-^H  coupling  constants  for  7  and  8  as  a  function  of  solvent  polarity.  16  The 
trend  of  increasing  coupling  constant  with  increasing  solvent  polarity  is  indicative  of  an  evolution 
towards  a  more  charge-separated  ground  state  (Figure  4,  right). 


RESULTS  AND  DISCUSSION 

The  solvent-dependent  polarizabilities  for  the  molecules  in  Figure  4  are  presented  in  Table  1. 
Molecules  1  and  2  display  positive  y  values  consistent  with  those  expected  given  the  large 
positive  bond  length  alternation  values  from  previous  structure  determinations.  The  increase  of 
y  with  increasing  conjugation  length  has  been  noted  by  other  authors.  Molecule  3,  with  the 
weak  formyl  acceptor,  displays  an  increasing  \L-^  product  and  a  positive  peak  in  y  with 
increasing  solvent  polarity.  This  behavior  is  consistent  with  that  expected  for  the  structural 
assignment  above  3  in  region  A  (Figure  2).  The  increase  in  y  for  the  donor-acceptor  polyene,  3, 
over  that  of  the  symmetrical  unsubstituted  polyenes  1  and  2  has  been  noted  previously^  ^  and  is 
due  to  the  lowering  of  symmetry  with  the  result  that  the  P2  term  (equation  3)  contributes 
positively  to  y. 

Compound  4  with  the  stronger  dicyanovinyl  moiety  displays  a  positive  peak  in  |i  |3  and  a 
sign  change  in  y.  In  fact,  the  peak  and  the  zero  crossing  of  y  occur  simultaneously, 
consistent  with  the  derivative-like  relationship  predicted  in  Figure  2  and  with  the  structural 
assignment  of  4  in  region  B. 

Molecules  5  and  6  display  decreasing  ^i•p  products  with  increasing  solvent  polarity.  For  5 
in  the  most  polar  solvent  and  for  6  in  moderate  polarity  solvents,  (i-P  changes  sign,  consistent 
with  the  ground-state  geometry  being  tuned  through  the  cyanine-limit  of  zero  bond  length 
alternation.  The  y  values  for  5  become  large  and  negative  with  increasing  solvent  polarity  with 
some  evidence  that  y  peaks  in  a  negative  sense  in  the  most  polar  solvents.  Molecule  6  displays  a 
clear  negative  y  peak  in  moderate  polarity  solvents.  The  zero  crossing  of  the  ^t  p  curve  for  both 
5  and  6  coincides  with  the  negative  y  peak,  consistent  with  the  proposed  structure-property 
relationships.  Additionally,  to  the  extent  to  which  a  can  be  represented  by  (l/Egg),  a  is 
maximized  at  the  zero  crossing  of  |i-P,  again  consistent  with  the  predicted  trend  (Figure  2). 

By  increasing  the  ground-state  polarization  further,  as  with  molecules  7  and  8,  there  is  some 
evidence  that  the  degree  of  mixing  required  to  maximize  |i-p  in  a  negative  sense  has  been 
achieved.  A  zero-crossing  in  y  for  7  in  nonpolar  solvents  is  observed  as  well  as  evidence  of  y 
approaching  a  positive  peak  for  8  in  the  most  polar  solvents.  The  zero  crossing  of  y  however 
does  not  coincide  with  the  negative  peak  in  |i-p.  This  may  be  explained  by  the  large  ground-state 
dipole  moment  shifting  the  peak  in  pi-p  to  lie  closer  to  the  y  peak  than  would  the  P  peak. 


Table  1.  Maximum  absorption  wavelength  (Xmax,  units  of  nm),  solvent-dependent  p.p  (units  of 
10-48  esu)  and  y  (units  of  10-36  esu)  for  1-8.  The  estimated  precision  in  p,-P  and  y  is  ±15%. 
The  polarity  of  the  solvents  increase  (the  contribution  of  the  charge-separated  structure  to  the 
ground-state  geometry  increases)  from  left  to  right.  Normalized  Ex(30)  values  16  of  the  solvents 


The  structure-property  relationships  proposed  in  Figure  2  can  be  further  tested  by  examining  the 
third-order  polarizability  of  molecules  possessing  zero  bond  length  alternation,  for  example 
cyanines  and  oxonols.  The  results  for  a  series  of  these  chromophores  are  presented  in  Table  11 
where  it  can  be  seen  that,  with  the  exception  of  the  shortest  chromophores,  the  y  values  are 

negative  as  predicted. 

Table  II.  y  values  (units  of  10‘36  esu)  and  maximum  absorption  wavelengths  (A-max’  of 
nm)  for  molecules  possessing  zero  bond  length  alternation.  The  estimated  precision  in  y  is  ±15%, 
except  in  the  case  where  the  magnitude  of  y  <  3  x  10'36  esu  when  the  estimated  error  in  y  is  ±5  x 
10"3^  esu. 


Molecule 

n 

1 

'^max 

Y 

cio7 

1 

318 

+  3 

2 

416 

-40 

3 

516 

-370 

CN 

1 

438 

-45 

NC  "+CN 

N(CH3)4 

2 

538 

-500 

/  \ 

1 

2 

266 

362 

+  3 

-7 

SUMMARY 


Donor-acceptor  polyenes  of  comparable  conjugation  length  have  been  synthesized  and  their 
solvent-dependent  hyperpolarizabilities  probed  by  EFISH  and  THG  measurements  at  1907  nm. 
Optimization  in  both  a  positive  and  negative  sense,  as  well  as  sign  changes  in  p-P  and  y  were 
observed.  These  observations  were  explained  by  molecular  structure  changes  (probed  by  x-ray 
crystallography,  1H-NMR,  electronic  absorption  and  Raman  spectroscopies)  resulting  from  the 
variation  of  mixing  of  neutral  and  charge-separated  resonance  forms  upon  changing  the 
donor/acceptor  strengths  and  solvent  polarity.  The  trend  of  the  geometry-dependent 
hyperpolarizabilities  is  fully  consistent  with  the  trends  predicted  using  geometry-optimized 
electric-field- dependent  AMI  and  INDO-CI  calculations  and  with  the  prediction  that  the 
polarizabilities  are  related  by  a  derivative  relationship.  A  separate  series  of  ionic  molecules 


possessing  zero  bond  length  alternation  were  found  to  exhibit  negative  third-order  polarizabilities, 
also  in  agreement  with  the  proposed  structure-property  relationships. 
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The  Effect  of  Molecular  Polarization  on  Bond  Length  Alternation,  Linear  Polarizahility, 
First  and  Second  Hyperpolarizability  in  Donor-Acceptor  Polyenes  as  a  Function  of 
Chain  Length 

Introduction 

Donor-acceptor  polyenes  are  the  prototypical  organic  chromophores  that  have  been 
considered  for  their  nonlinear  optical  (NLO)  response.  In  attempting  to  understand  the 
relationship  between  chemical  structure  and  the  linear  and  nonlinear  optical  properties  of 
molecules,  two  relevant  molecular  parameters  are  typically  varied;  ground-state  polarization  (by 
varying  the  relative  strengths  of  the  donor  and  acceptor),  and  the  molecular  length.  We  have 
recently  described  a  procedure  that  permits  us  to  calculate  a,  P  and  y  of  donor-acceptor  polyenes 
of  the  form  R2N-(CH=CH)n-CHO  ,(1)  through,  a'range  of Tond"  length  alternation  (BLA) 
corresponding  to  differing  degrees  of  charge-separation  in  the  molecule  (Figure  2,  shown  here  for 
n=5  and  previously  reported  for  n=A)}‘  ^  Using  this  method,  the  molecule  can  be  varied  from  a 
gas-phase  form  with  little  charge  separation  [the  polyene  limit]  depicted  by  the  resonance  form  on 
the  left  side  of  Figure  1,  through  a  bond  length  equivalent  form  in  which  both  resonance  forms 
make  an  approximately  equal  contribution  [the  cyanine  limit],  through  a  strongly  charge-separated 
form  [the  zwitterionic  limit]  dominated  by  the  resonance  structure  at  the  right  side  of  Figure  1. 

Having  established  the  basic  validity  of  this  approach,  we  now  present  results  of  a  computational 
study  that  varies  both  BLA  and  molecular  length.  We  point  out  that  a  physically  accessible 
parameter,  BLA,  can  be  related  to  maxima  and  minima  in  the  molecule's  linear  polarizability  (a), 
and  first-  (p)  and  second  molecular  hyperpolarizability  (y)  for  molecules  of  differing  length.  In 
addition,  we  provide  a  more  detailed  description  of  the  computational  technique  as  well  as  several 
tests  to  help  validate  the  procedure. 

Results  and  Discussion 

Using  the  semi-empirical  AMI  parameterization  in  the  MOP  AC  package,^' we  examined 
donor-acceptor  polyenes  of  the  form  1  (R  =  Me,  n  =  2-7)  under  the  influence  of  an  external 
perturbation  designed  to  vary  the  ground-state  polarization  and  geometry.  This  perturbation  was 
two  positive  and  two  negative  point  charges  placed  as  is  shown  in  Figure  1  and  moved 
incrementally  towards  the  molecules  (from  40  A  to  4  A).  While  the  distances  between  the  point 
charges  and  the  molecule  were  held  fixed,  the  molecular  geometry  was  optimized.  As  the  point 
charges  were  placed  closer  to  the  molecules,  it  was  observed  that  the  single  bonds  in  the 
molecules  became  shorter  (from  =  1.44  A  to  -  1.35  A)  and  the  double  bonds  became  longer 
(from  »  1.35  A  to  =  1.44  A).  Although  this  entire  range  of  bond  lengths  could  not  be  spanned 


Figure  1.  The  neutral  (top,  left)  and  charge-separated  (top,  right)  resonance  structures  that 
contribute  to  both  the  electronic  and  nuclear  structure  of  the  molecule  and  a  depiction  of  the 
placement  of  the  point  charges  (bottom).  In  the  absence  of  an  external  perturbation,  a 
computational  geometry  optimization  is  biased  against  charge  separation  and  depicts  the  molecule 
almost  solely  as  it  is  represented  on  the  left.  Point  charges  promote  increased  charge-separated 
character  as  they  approach  the  molecule. 

for  the  shorter  polyenes  under  study  (n  <  4),  the  range  corresponds  to  a  change  in  bond  length 
alternation  (defined  as  the  difference  between  the  average  of  the  central  n-2  C=C  bonds  and  the 
average  of  the  central  n-2  C-C  bonds  in  the  molecules  with  the  resonance  structure  shown  at  the 
left  of  Figure  1)  from  BLA  =  -0.09  A,  represented  predominantly  by  the  charge-neutral  (left) 
resonance  structure  in  Figure  1  to  BLA  =  +0.09  A,  represented  predominanUy  by  the  charge- 
separated  (right)  resonance  structure  in  Figure  1.  At  the  same  time,  the  ground-state  dipole 

moments  of  the  molecules  increased  smoothly. 

Values  for  a,  P  and  y  for  the  polarized  molecules  were  then  calculated  using  a  finite  field 

subroutine.^  Although  polarization  by  point  charges  is  different  than  charge  stabilization  by 
solvent  (in  particular,  the  largest  fields  applied  in  our  calculations  are  considerably  larger  than  the 
reaction  fields  for  common  solvents),  this  method  qualitatively  reproduces  experimental  trends  in 
geometry  and  polarizabilities  as  a  function  of  increasing  ground-state  polarization.^' Several 
observations  are  made  from  Figure  2.  (i)  The  value  of  a  exhibits  a  peak  at  the  cyanine  limit,  (ii) 
The  value  of  P  exhibits  a  positive  peak  between  the  polyene  and  the  cyanine  limits,  crosses 
through  zero  at  roughly  the  cyanine  limit  (at  which  point  the  difference  between  the  dipole 
moment  of  the  strongly  allowed  first  excited  state  and  the  ground  state  is  roughly  zero22),  and 
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exhibits  a  negative  peak  between  the  cyanine  and  zwitterionic  limits,  (iii)  Between  the  polyene 
and  cyanine  Hmits,  the  value  of  y  is  initially  positive  and  increases,  reaching  a  peak.  The  value  of 
ythen  decreases,  crosses  through  zero  (at  a  BLA  corresponding  to  that  where  P  peaks),  becomes 
negative,  reaching  a  negative  peak  (with  a  maximized  absolute  value)  at  the  cyanine  limit  To  first 
approximation,  the  region  between  the  cyanine  and  zwitterionic  limits  mirrors  that  on  the  other 
half  of  the  curve,  (iv)  The  curves  are  roughly  derivatives  with  respect  to  each  other  as  is 
discussed  in  detail  elsewhere.^^  For  P  and  y,  solvent  dependent  electric  field-induced  second 
harmonic  generation  (EFISH)  and  third  harmonic  generation  (THG)  studies  on  a  series  of 
molecules  of  similar  length  and  with  ground-state  polarizations  that  span  almost  the  entire  range  of 
BLA  depicted  by  these  curves  confirm  the  behavior  predicted  in  Figure  2.®'  Furthermore,  it  has 
recently  been  shown  that  molecules  with  reduced  BLA  similar  to  that  at  the  pocti'-’c  in  the  3 
curve  have  exceptionally  high  nonlinearities.^^  As  discussed  earlier,’'  ‘  ail  of  these  results  am 
consistent  with  simplified  perturbative  expressions  which  approximate  a,  p,  and  y  based  upon 
couplings,  transition  energies,  and  changes  in  dipole  moments  between  the  ground  and  one  or 
two  excited  states  of  the  molecule.^' In  addition,  recent,  similar  computational  studies  using 
a  sum-on-states  approach,^^  using  an  anharmonic  oscillator  model,^^  and  using  a  simplified  two 
level  valence  bond  approach^^  predict  essentially  the  same  relationships  as  a  function  of  ground- 
state  polarization. 
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Figure  2.  Correlation  of  molecular  polarizability,  a  (top  plot),  first ,  P  (center  plot)  and  second, 
Y  (bottom  plot )  hyperpolarizabilities  (on  the  left  axes),  with  bond  length  alternation  (BLA, 
defined  in  the  text)  for  Me2N-(CH=CH)5-CHO  (•).  The  bond  length  alternation  at  which  a,  (top 
plot),  P  (center  plot)  and  y  (bottom  plot)  are  maximized  (♦),  are  minimized  (■)  or  cross  through 
zero  (A)  for  a  given  of  chain-length  (n),  where  n  is  shown  on  the  right  axes.  The  shorter 
molecules  could  not  be  driven  to  the  most  charge  separated  form  (BLA  >  =  4-0.025),  so  not  all 
maxima/minima  could  be  reported.  Likewise,  the  7-double  bond  polyene  could  not  be  driven 
smoothly  through  the  region  of  BLA  =  0. 
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When  evaluating  the  molecular  length  dependence  of  j3,  it  is  common  to  fit  data  to  a 
function  of  the  form  p(0)  N“  25-28  Here,  N  is  defined  as  the  number  of  conjugated  multiple 
bonds  in  the  molecule  (including  the  C=0  bond,  i.e.  N  =  n  +  1)  so  as  to  more  accurately  represent 
the  real  length  of  the  molecule.  The  exponent,  m,  lends  insight  into  the  chain  length  dependence 
of  p.  We  have  computed  m  for  P  for  the  molecules  at  various  bond  length  alternations  and  for 
molecules  at  Pmax  (Table  1).  The  exponent  m  is  also  computed  for  a  and  y.  Exponents  derived 
from  data  in  the  literature  are  also  reported  (for  these  cases  N  is  defined  as  above  and  therefore 
exponents  may  differ  from  those  reported  in  the  original  work,  due  to  the  use  of  altnerative 
counting  scheme  for  the  conjugation  length).^^  The  exact  bond  length  alternation  corresponding 
to  a  maximized  P  and  y  varies  somewhat  with  molecular  length,  and  this  value  is  graphically 
displayed  as  a  function  of  chain  length,  n,  in  Figure  2._  It  is  not  expected  that  each  member  of  this 
homologous  series  of  molecules  will  have  the  same  bond  length  alternation  in  a  given  solvent,  but 
the  range  of  exponents,  for  both  constant  bond  Idrr^th'kltemation  and  for  bond  length  alternation 
corresponding  to  a  maximized  p  or  y  gives  a  range  into  which  optimized  molecules  are  expected  to 
fall.  In  particular,  it  is  expected  that  increasingly  strong  donors  and  acceptors  will  be  required  to 
reach  a  given  BLA  value  as  molecules  become  longer  reflecting  the  greater  difficulty  in  separating 
charge  over  greater  lengths.  The  exponents  at  or  near  the  maxima  or  minima  of  the  as  ,  P,  and  y, 
curves  are  larger  than  those  computed  for  hyperpolarizabilities  based  upon  bond  length  alternate 
(gas-phase)  geometries,^^'  28  and,  in  the  case  of  p,  for  those  observed  for  molecules  containing 
aromatic  endgroups.^'  50-32  por  a,  the  data  shows  that  chromophores  with  minimal  bond  length 
alternation  (the  cyanines)  are  optimized,  a  conclusion  consistent  with  previous  experimental 
determinations.55  p  is  optimized  at  BLA  =  -0.05±0.01  A,  and  around  this  point,  m  «  4.5.  We 
have  synthesized  simple  donor-acceptor  polyenes  such  as  1  (n  =  0-3)  as  well  as  R2N-(CH=CH)n- 
CHC(CN)2,  2  (n=0-3)  and  find  that  the  value  for  m  for  these  molecules  (measured  by  EFISH  in 
CHCI3)  is  3.2  and  5.4,  respectively.^^  As  mentioned  earlier,  a  series  of  vinylogous  compounds 
(such  as  2)  will  not  necessarily  have  the  same  BLA  at  different  chain  lengths.  Therefore,  it  is 
possible  to  move  closer  to  the  maxima  of  the  curve  with  increasing  chain  length,  if  BLA  for  a 
short  compound  in  the  series  is  too  close  to  zero  to  optimize  p  and  if  IBLAI  becomes  larger  with 
increasing  length.  This  would  result  in  higher  exponents  than  would  be  predicted  for  any  given 
bond  length  alternation.  Conversely,  if  the  magnitude  of  BLA  is  greater  than  that  required  to 
optimize  P  in  a  positive  sense,  then  for  a  given  donor  or  acceptor,  the  exponent  is  expected  to  be 
lower  than  might  be  predicted  based  on  short  compounds  within  the  series.  The  magnitude  of  y  is 
maximized  in  the  cyanine  limit,  but  also  shows  a  maximum  at  BLA  =  -0.07±0.01  A  which  may 
correspond  to  more  stable  molecules,  particularly  for  longer  chain  lengths.  The  exponents 


corresponding  to  these  extrema  are  again  significantly  greater  than  those  predicted  previously  for 
polyenes.^^' 

Table  1.  Summary  of  exponents  for  a,  p,  and  y  of  polyenes  and  polarized  polyenes _ 


Quantity 

Geometry 

m 

N's  Used 

Compounds/Methodi’ 

a 

BLA  =  -0.09 

1.33 

3-8 

This  work 

Maximum* 

2.16 

3-8 

This  work 

Computational 

1.38 

3-8 

H-(CH=CH)n-H.  N  =  n  using  the 

Reference  ^ 

6-3 IG  +  PD  basis  set  w.  SOS 

Computational 

1.31 

3-8 

H-(CH=CH)e-H,  M  =  n  using 

Reference 

AMl/FF 

p 

BLA  -  -0.09 

2.95 

3-8 

This  work 

Maximum* 

4.51 

3-8 

This  work 

BLA  =  -0.06 

4.28 

3-8 

This  work 

BLA  =  -0.05 

4.49 

3-8 

This  work 

BLA  =  -0.04 

4.61 

3-8 

This  work 

Computational 

3.43 

3-7 

Me2N-(CH=CH)n-N02 

Reference 

N  =  n  +  1  using  CNDOVSB/SOS 

Computational 

3.2  F 

3,5,7,9 

Me2N-(CH=CH)n-CHO 

Reference 

N  =  n  +  1  using  CNDOVSB/SOS 

Computational 

3.52 

3-6 

Me2N-(CH=CH)n-N02, 

Reference 

N  =  n  +  1  using  PPP/SOS 

Experimental 

4.02 

4-6 

Me2N-Ph-(CH=CH)n- 

Reference^^ 

CH=C(CN)2,  N=n+4,  EFISH  at 

1.06  |im  in  chloroform. 

Experimental 

2.27 

6-9 

Me2N-Ph-(CH=CH)n-Ph-N02, 

Reference^^ 

N=n-i-5,  EFISH  at  1.06  |im  in 

chloroform. 
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Experimental 

2.48 

5,8,10,12 

^iP(O)  values,  Julolidine- 

Reference^^ 

(CH=CH)n-CH=C(a^)2,  N=a+4, 
EFISH  at  1.91  in  chloroform  or 

acetone. 

Experimental 

3.26 

2,5,7,9 

11(3(01  values,  Benzodithia- 

Reference^^ 

(CH=CH)n-CHO,  N=n+1,  EFISH 
at  1.34  |im  in  chloroform  or 

acetone. 

Experimental 

3.2 

1-4 

Me2N-(CH=CH)n-CHQr  N-n  1 1  ’ 

RefereiiCf^" 

EFISH  Ri  i  |ini  in  chMi^'foim 

Experimental 

5.4 

2-5 

Me2N-(C::-CI  vn  -CH=C(CN)2, 

Reference^ 

N=n+2,  EFISH  at  1.91  p,m  in 

chloroform 

Y  BLA  =-0.09 

4.400 

3-8 

This  work 

Maximum^ 

6.091 

3-8 

This  work 

BLA  =  -0.08 

5.111 

3-8 

This  work 

BLA  =  -0.07 

5.969 

3-8 

This  work 

BLA  =  -0.06 

6.727 

3-8 

This  work 

Minimum^<J 

6.509 

3-7 

This  work 

Computational 

3.421 

3-8 

H-(CH=CH)n-H,  N  =  n  using  the 

Reference  ^ 

6-3 IG  +  PD  basis  set  w.  SOS 

Computational 

4.004 

3-8 

H-(CH=CH)n-H,  N  =  n  using 

Reference 

AMl/FF 

aThat  bond  length  alternation  that  maximizes/rainiraizes  the  quantity  of  interest.  In  the 
case  of  Y,  the  maxima  refers  to  that  for  the  least  charge  separated  form  (i.e.  BLA  <  0). 
bComputational  results  are  for  co  =  0  (Sum-over-states  [SOS])  or  static  perturbing  field  (Finite 
Field  [FF]).  Experimental  results  for  (3  are  by  EFISH  =  electric  field-induced  second  harmonic 
generation.  Definition  of  N  is  consistent  throughout  and  is  different  from  that  used  by  other 
authors.  Two  double  bonds  are  counted  for  a  benzene  ring  and  an  additional  double  bond  is 
counted  for  the  bond  in  the  acceptor,  cpor  a  curve  containing  only  3  points  (N  =  3,5,7), 
m  =  3.281.  dCurve  fit  is  to  lYminl- 
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In  conclusion,  these  results  provide  structure-property  relationships  for  donor-acceptor 
polyenes,  both  as  a  function  of  molecular  length  and  as  a  function  of  the  ground-state  polarization 
and  resulting  bond  length  alternation  in  the  molecules.  In  particular,  we  predict  that  molecules 
with  a  significantly  lower  magnitude  bond  length  alternation  than  is  observed  in  polyenes  and 
diphenyl-capped  polyenes  can  have  larger  hyperpolarizabilities  as  well  as  faster  increases  in 
hyperpolarizability  with  increasing  chain  length.  Thus,  we  suggest  that  cyanines  with  0  A  of 
bond  length  alternation  and  donor-acceptor  polymethines  with  roughly  10.071  A  of  bond  length 
alternation  be  considered  for  third-order  NLO  applications.  Likewise,  we  suggest  that  donor- 
acceptor  polymethines  with  roughly  10.051  A  of  bond  length  alternation  be  considered  for  second- 
order  electrooptic  applications. 

Description  of  the  computation. 

General 

Geometry  optimizations  of  1  (R  =  Me,  n  =  2-7)  were  accomplished  using  the  AMI 
parameterization  available  in  MOP  AC  6.0.^'  ^  Using  a  Z-matrix  definition  for  the  starting 
molecular  geometries,  the  distances,  angles  and  dihedral  angles  between  external  point  charges  of 
opposite  signs,  and  between  two  of  the  point  charges  (positive  Sparkles)  and  the  molecule 
(nitrogen  atom)  were  fixed.  To  accomplish  this  task  the  molecules  were  held  planar  during  this 
procedure.  Using  a  less  convenient  external  point  charge  arrangement  in  which  the  point  charges 
were  placed  in  the  molecular  plane,  it  was  shown  that  p,,  a,  P,  and  y  versus  bond  length 
alternation  were  basically  unchanged  when  the  least  planar  molecule,  1  (n  =  3),  was  alternately 
held  planar  and  permitted  to  become  nonplanar  during  geometry  optimization.  Likewise  these  two 
arrangements  of  the  external  point  charges  around  the  molecule  (e.g.  point  charges  in  the 
molecular  plane  versus  point  charges  perpendicular  to  the  molecular  plane)  made  little  difference 
in  p,,  a,  p,  and  ywhen  these  quantities  were  plotted  versus  bond  length  alternation.  The  detailed 

shape  of  the  electric  field  will  of  course  be  dependent  upon  the  arrangement  of  the  external  point 
charges  (Figure  1).  This  point  has  been  considered,  and  alternate  arrangements  of  the  point 
charges  were  employed.  These  included  a)  four  charges  in  the  molecular  plane,  b)  eight  charges 
arranged  as  dipoles  perpendicular  to  the  molecular  plane,  and  c)  four  charges  at  different  distances 
from  the  donor  and  acceptor  to  account  for  donors  and  acceptors  of  different  relative  strengths. 
For  all  of  these  different  arrangements,  and  even  when  the  potential  stabilization  at  one  side  of  the 
molecule  was  3  times  that  of  the  other  (corresponding  to  a  much  weaker  donor  than  acceptor  and 
vice  versa)  the  behavior  of  p,  a,  p,  and  y  were  essentially  unchanged,  and,  in  most  cases,  curves 
were  superimposable  when  viewed  versus  bond  length  alternation.  Again,  it  is  not  our  intention 
to  simulate  the  effect  of  a  particular  solvent  or  other  medium  upon  a  molecule  with  any  given 
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external  point  charge  arrangement  but  rather  to  probe  the  electronic  and  geometric  changes  in  the 
molecule,  and  its  (hyper)polarizabilities  as  it  becomes  more  polarized.  In  this  case,  these  results 
merely  indicate  that  the  specific  point  charge  arrangement  used  to  polarize  the  molecule  is  not 
critical  to  observing  these  effects. 

The  point  charges  (Sparkles)  were  moved  in  steps  from  40  A  to  4A  from  the  ends  of  the 
molecule  in  the  trajectory  shown  in  Figure  1.  The  step  increment  became  smaller  as  the  point 
charges  approached  the  molecule  and  the  change  in  electric  field  with  a  change  in  distance  became 
larger,  beginning  with  an  increment  of  10  A  (40  A,  30  A,  20A),  followed  by  an  increment  of  1  A 
(15  A  through  8  A),  and  ending  with  an  increment  of  0.2  A  (8  A  through  4  A)  for  n  =  3-6.  For  n 
=  7,  this  last  increment  was  varied  as  little  as  0.0002  A  (between  6.8820  A  and  6.8840  A)  in  ah 
effort  to  map  out  the  region  around  BLA  =  0  A.  Despite  the  use  of  this  small  increment,  this 
region  was  not  mapped  out.  BLA  jumped  from  -0.013  A  (for  a  Sparkle  distance  of  6.8840A)  to 
+0.015  A  (for  a  Sparkle  distance  of  6.8838  A).  Reasons  for  the  existence  of  this  "forbidden 
region"  in  molecules  of  this  length  (and  longer,  results  not  presented  here)  are  under  consideration 
but  may  be  related  to  charge  localization  in  infinite  polyenes^^'  and  long  cyanine 
molecules.^^' 

At  each  fixed  Sparkle  distance,  the  geometry  was  optimized,  and  |i,  a,  P  and  y  were 
calculated  in  the  presence  of  these  point  charges  using  the  finite-field  subroutine^  available  in 
MOP  AC  6.0  using  an  SCF  convergence  criterion  (SCFCRT)  of  lO'^O  a.u.  as  recommended  by 
Kurtz.^  All  of  these  quantities  are  vectors  or  tensors,  and,  in  most  cases,  average  values  are 
reported.  These  average  values  are:  a  =  (ttxx  +  ocyy  cczz)/3;  p|i,  =  (3/5)*(P*p/llil)  where  P»p.  = 
Pxl^x  +  PyM-y  +  Pzl^z  and  p/  =  {P,//  +  ^ijj+  P/jt/:};  and  y  =  (1/5)*[Yxxxx  +  Yyyyy  +  Tzzzz  + 
2.0*(yxxyy  +  Yxxzz  +  Yyyzz)]-  All  reported  values  are  obtained  from  the  energy  equations^  (E4) 
in  the  finite  field  subroutine.  Values  obtained  from  the  dipole  equations^  (DIP)  differ  by  less  than 
1%  of  the  E4  values. 

It  is  important  to  note  that  although  two  explicit  fields  are  applied  to  the  molecule,  the 
point  charge  field  and  the  electric  field  within  the  finite  field  calculation,  these  fields  have 
distinctly  different  effects  on  the  molecule.  In  particular,  the  first  field  is  permitted  to  influence 
both  its  geometric  and  electronic  structure.  Thus,  at  the  start  of  the  finite  field  computation,  the 
molecule  is  in  a  new  equilibrium  geometry  and  electronic  configuration,  within  which  the 
hyperpolarizabilities  are  then  calculated  (Figure  3).  By  varying  the  strength  of  the  electric  field 
employed  in  the  finite  field  routine  both  numerical  accuracy  problems  and  electron  configurational 
changes,  observed  for  field  strengths  that  are  two  small  and  large,  respectively  were  avoided  (see 
below).  Using  point  charges  to  polarize  the  electrons  without  permitting  the  molecular  geometry 
to  change  results  in  a  finite  field  calculation  of  a  molecule  in  a  non-equilibrium,  polarized,  state. 
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Likewise,  removal  of  the  point  charges  after  geometry  optimization,  but  before  the  finite  field 
procedure  leaves  nuclei  in  a  position  corresponding  to  a  polarized  molecule,  but  removes  the 
potential  that  stabilizes  the  polarization.  Both  of  the  above  procedures  were  performed;  however, 
for  each  case,  since  the  nuclei  and  the  electrons  were  not  in  equilibrium  configurations,  much 
larger  applied  fields  were  needed  to  polarize  the  molecule,  as  compared  to  applying  a  field  and 
allowing  the  geometry  to  optimize.  These  results  suggested  to  us  that  both  the  applied  electric 
field  and  the  position  of  the  nuclei  make  important  contributions  to  the  polarization  response  and 
that  the  most  reasonable  procedure  was  to  account  for  both  effects. 


24 


o 


Figure  3.  Molecule  1  [top]  in  the  absence  of  an  external  electric  field,  [middle]  exhibiting  both 
nuclear  and  electronic  polarization  due  to  external  point  charges,  and  [bottom]  exhibiting  further 
electronic  polarization  (indicated  by  the  arrows)  as  the  result  of  the  finite  field  procedure. 

Tests  of  the  finite  field  procedure 

In  the  finite  field  procedure,  the  default  strength  of  the  perturbing  electric  field  (10*3  a.u.) 
was  employed,  and  all  values  reported  are  those  obtained  using  this  base  field  strength.  To  test 
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the  validity  of  this  field  strength,  it  was  varied  to  ensure  that  both  numerical  accuracy  problems 
and  electron  configurational  changes,  observed  for  field  strengths  that  are  too  small  and  large, 
respectively  were  avoided.®  Figure  4  shows  the  value  of  y  (=  (1/5)*[Yxxxx  +  Yyyyy  ■*"  Yzzzz  + 
2*(Yxxyy  +  Yxxzz  +  Yyyzz)])  for  1  (n  =  4)  both  in  the  presence  and  absence  of  the  external  point 
charges  arranged  as  in  Figure  1.  This  longer,  more  (hyper)polarizable  molecule  is  much  more 
sensitive  to  changes  in  applied  field  than  were  the  smaller  substituted  benzenes  reported  earlier  by 
Kurtz,®  indicating  that  this  test  is  worth  performing  for  these  molecules.  However,  the  external 
point  charges  do  not  significantly  change  the  sensitivity  of  the  computed  (hyper)polarizabilities  to 
the  applied  field.  This  observation  is  consistent  with  the  conclusion  that  the  molecule  is  in  an 
equilibrium  nuclear  and  electronic  configuration  at  the  start  of  the  finite  field  procedure. 


•j  000  H - 1 - ' — — I  I  I'M - ' - ' — ' — 1  I  'll 


-1  oooH - ^ — ' — '  '  ‘  ■ ‘-H - ^ — ' — '  '  ‘  '  I 

10"^  0.001  0.01 

Applied  field  (a.u.) 


Figure  4.  The  sensitivity  of  y  (report  ed  Figure  2),  calculated  using  energy  equations®  (E4)  values 
from  the  finite  field  subroutine,  for  1  with  [circles]  no  point  charges,  [triangles]  external  point 
charges  placed  10  A  from  the  molecule  as  in  Figure  1  and  [boxes]  external  point  charges  placed  6 
A  from  the  molecule  as  in  Figure  1.  Over  the  range  shown,  the  dipole  equation  values  were 
comparable  to  those  shown  down  to  a  field  strength  of  2  x  lO-^  a.u.  at  which  point  they  diverged. 
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Results  at  non-equilibrium  configurations 

Figure  5  shows  p  computed  for  2  in  what  we  rationalize  as  non-equilibrium 
configurations.  In  one  case,  (Figure  5,  top)  the  molecule  was  geometry  optimized  in  the  gas 
phase  and  external  point  charges  were  placed  around  the  molecule  only  for  the  finite  field 
calculations,  corresponding  to  a  computation  in  a  non-equilibrium  configuration.  The  value  of  P 
increases  as  the  molecule  is  polarized  and  even  maximizes  at  approximately  the  same  induced 
dipole  moment  as  that  observed  when  the  geometry  is  permitted  to  optimize.  The  absence  of 
nuclear  relaxation,  however,  limits  the  extent  of  molecular  polarization.  In  the  second  case, 
(Figure  5,  bottom)  the  molecule  was  geometry  optimized  in  the  presence  of  the  external  point 
charges,  and  these  were  removed  before  the  finite  field  calculation,  corresponding  to  another  non¬ 
equilibrium  configuration.  The  molecule  polarizes,  as  evidenced  by  an  increase  in  dipole  moment 
from  6  debye  to  >  30  debye  over  the  range  shown  in  the  graph.  This  polarization  again  makes  an 
incomplete  contribution  to  the  overall  molecular  polarization,  evidenced  by  an  even  more  sluggish 
response  to  the  electric  field  applied  in  the  first  half  of  the  procedure.  Overall,  it  is  our  contention 
that  both  electronic  and  nuclear  polarization  make  an  important  contribution  to  the  configuration  of 
the  molecule,  and  both  should  be  taken  into  account. 


Figure  5.  P  calculated  in  two  non-equilibrium  scenarios  for  2.  [Top]  Plot  of  P  versus  jJ.  where 
the  molecule  was  geometry  optimized  in  the  gas  phase  and  external  point  charges  were  placed 
around  the  molecule  (as  in  Figure  1)  only  for  the  finite  field  calculations,  corresponding  to  a 

computation  in  a  non-equilibrium  nuclear  configuration.  [Bottom]  Plot  of  P  versus  bond  length 
alternation  the  molecule  was  geometry  optimized  in  the  presence  of  the  external  point  charg^,  and 
these  were  removed  before  the  finite  field  calculation,  corresponding  to  a  non-equilibrium 

electronic  configuration.  P  is  plotted  versus  bond  length  alternation. 
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The  Effects  of  Aromaticity  and  Reduced  Bond  Length  Alternation  on  Linear  and 
Third-Order  Polarizabilities  of  Conjugated  Organic  Chromophores 

Introduction 

The  induced  polarization  of  a  molecule  by  an  applied  electric  field,  E,  can  be  approximated 
by  the  power  series  expansion: 

Pinduced  =  ttijE  +  Bijk/2  E*E  +  Yykl/6E»E*E  +  ...  (1) 

Relationships  between  molecular  structure  and  the  linear  polarizability  (a),  first  hyperpolarizability 
(p)  and  second  hyperpolarizability  (y)  are  of  interest  in  the  study  of  organic  materials  for  nonlinear 
optics  (NLO).  In  this  paper,  we  employ  a  computational  method  (AMI  coupled  with  finite  field)  to 
examine  these  relationships  for  two  classes  of  organic  molecules:  bipolaronic  model  compounds 
(organic  dications/dianions)  and  cyanines  (Figure  1).  We  rationalize  the  results  of  these 
computations  by  considering  the  effects  of  aromaticity,  charge  localization  and  bond  length 
alternation  upon  the  various  molecules. 

Procedure 

Using  the  semi-empirical  AMI  parameterization  in  the  MOP AC  package,^’  ^  structures  1-5 
(Figure  1)  were  geometry-optimized  in  the  neutral  (uncharged)  form,  in  the  dication  form  for  1, 3 
and  4  and  in  the  dianion  form  for  2  and  5.  Cationic  cyanine  6  was  also  geometry-  optimized  for 
different  molecular  lengths  (n  =  2-10).  The  geometry-optimization  was  not  conducted  with  any 
defined  symmetries,  and  for  1-5,  molecular  geometries  were  found  to  be  mirror  symmetric 
perpendicular  to  the  long  axis  of  the  molecule  within  0.01  A.  After  geometry  optimization,  a  and  y 
were  calculated  for  all  by  a  finite  field  subroutine.^  Average  values  for  these  tensors  are  reported, 
with  a  =  (l/3)*(axx  +  “yy  +  “zz)  and  y  =  (l/5)*[yxxxx  +  Yyyyy  +  Yzzzz  +  2.0«(yxxyy  +  Yxxzz  + 
yyyzz)]-  Values  reported  (Table  I)  are  those  computed  using  the  energy  (E4)  equations  from 
MOP  AC.  Those  computed  using  the  dipole  (DIP)  equations  differ  by  less  than  2  %  from  the  E4 
values. 
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Figure  L  Molecules  investigated  in  this  study. 
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Table  I.  Linear  (a)  and  third-order  (7)  polarizabilities  for  1-6. 


1  10.53  782.3  2  5.29  -124.9 

1  (2+)  17.13  2595.0  3  8.51  -657.2 

2  (2.)  11.89  1633.3  4  13.03  -2439.3 

2  10.85  1070.9  5  18.87  -7207.0 

3  (2+)  10.15  1120.5  6  26.41  -18295.9 

3  11.69  1264.3  7  32.39  -14652.9 

4  11.59  1354.3  8^  28.34  30758.7 

4  (2-h)  15.32  3137.1  9^  29.25  47736.9 

5  8.55  726.6  10^  28.65  44251.7 

5  (2-)  14.84  1938.9 

a  X  10-23  esu;  ^  x  10-36  esu;  ^Charge  localization/bond  length 

alternation  observed  —  see  Figure  3. 

EFFECT  OF  OXIDATION/REDUCTION  (BIPOLARON  FORMATION)  ON  a  and 
Y  OF  CONJUGATED  ORGANICS 

Increased  y  has  been  predicted  upon  bipolaron  formation  (oxidation  or  reduction  by  two 
charges)  for  a  number  of  conjugated  organic  systems,  including  polyene-derived  chromophores.'*’ 
5  We  explore  this  hypothesis  for  compounds  1-5,  and  attempt  to  correlate  changes  in  a  and  y  with 
concomitant  changes  in  bond  length  alternation  and  aromaticity  upon  oxidation  or  reduction. 
Model  compounds  1-5  are  intended  to  represent  synthetically  accessible,  discrete  molecules  that 
have  different  molecular  topologies  and  that  are  sufficiently  long  to  allow  for  significant  charge 
separation.^’  ^ 

In  our  computations,  polyene-based  compounds  4  (upon  oxidation)  and  5  (upon 
reduction)  display  the  predicted  increases  in  a  and  y.  As  expected,  in  the  bipolarons,  the  charges 
build  up  near  the  ends  of  the  molecules  due  to  their  mutual  coulombic  repulsion.  Concurrently, 
bond  length  alternation  (defined  as  the  difference  between  the  lengths  of  the  single  and  double 
bonds,  where  only  non-ring  bonds  are  considered,  vide  infra)  is  reduced  near  the  ends  of  the 
molecule  (Table  II),  presumably  due  to  partial  delocalization  of  the  charges.  In  longer  systems, 
particularly  conjugated  polymers,  further  reduction  in  bond  length  alternation  due  to  charge 
delocalization  might  be  observed  since  the  mutual  coulombic  repulsion  may  be  greatly  diminished 
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as  the  charges  can  separate  and  therefore  delocalize  over  a  greater  distance.  Properties  such  as  the 
linear  absorption  spectra,  solvatochromism,  and  linear  and  hyperpolarizabilities  of  organic 
chromophores  are  sensitive  to  changes  in  bond  length  alternation  associated  with  the  introduction 
of  formal  charges  or  increased  charge  separation.*'^ ^  Such  effects  have  been  discussed  when 
comparing  the  optical  spectra  of  polyenes,  R-(CH=CH)n-R,  merocyanines,  R2N-(CH=CH)n.i- 
CH=0  and  cyanines,  R2N-(CH=CH)n.i-CH=NRf  of  similar  molecular  lengths.*  Also,  mid-gap 
(carrier)  states  in  conjugated  polymers  formed  upon  doping  (oxidation  or  reduction)  alter  their 
geometry,  changing  their  electronic,  and  linear  and  non-linear  optical  properties. Recently, 
enhanced  y  has  been  observed  in  a,0)-diphenyl-octatetraene  upon  photoexcitation  using  a  pump- 
probe  experiment. Although  the  equilibrium  geometry  of  this  excited-state  has  not  been 
determined,  molecular  orbital  descriptions  of  related  polyenes  invoke  reduced  bond  order 
alternation  in  the  first  allowed  excited-state.''*’ 


Table  II.  Selected  bond  lengths  (r,  in  A,  according  to  labelling  scheme  below)  for  the 
neutral/charged  forms  of  4  and  5 


Y 

Cmpd 

a 

b 

c 

d 

e 

f 

g . 

h 

(CH3)2N 

4 

1.39 

1.36 

1.44 

1.35 

1.44 

1.35 

1.44 

1.35 

(CH3)2N+ 

4(2+) 

1.33 

1.43 

1.37 

1.43  . 

1.36 

1.44 

1.35 

1.45 

0 

5 

1.23 

1.47 

1.34 

1.45 

1.35 

1.44 

1.35 

1.44 

o- 

5(2-) 

1.26 

1.41 

1.39 

1.38 

1.42 

1.36 

1.43 

1.35 

In  addition  to  effects  arising  from  changes  in  bond  length  alternation,  we  find  that  the 
topology  of  a  system  containing  aromatic/quinoidal  rings  is  important.  For  example,  when 
aromaticity  increases  upon  oxidation  (as  with  3  =>  3  (2+) ),  y  actually  decreases.  As  indicated  in 
Figure  1,  1,  2,  and  3  were  considered  both  in  neutral  and  in  charged  forms.  Resonance 
structures  shown  in  Figure  1  approximately  reflect  the  observed  geometry-optimized  structures. 
Selected  bond  lengths  are  shown  in  Table  IE. 

In  1,  the  molecule  is  aromatic  in  the  uncharged  state  and  becomes  quinoidal  upon 
oxidation.  Its  bond  lengths  suggest  that  the  charges  are  localized  near  the  ends  of  the  molecule. 
The  N-Csp2  bond  length  decreases  from  1.41  A  to  1.35  A  upon  oxidation,  consistent  with 
increased  bond  order  and  assumption  of  positive  charge  by  the  nitrogen.  Reduced  bond  length 
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alternation  is  observed  immediately  adjacent  to  the  ring  (bonds  d  and  e.  Table  HI)  compared  to  that 
observed  in  the  center  of  the  chain  (r(C=C)  =  1.35  A  and  r(C-C)  =  1.44  A  for  the  central  double 
and  single  bonds  in  1(2+),  respectively).  Reduced  aromaticity  in  the  oxidized  form  is  also 
expected  to  increase  a  and  y.  In  the  aromatic  form,  the  ^-electrons  are  effectively  delocalized 
within  the  rings;  however,  poor  coupling  between  the  rings  and  the  adjacent  polyene  prevent 
effective  delocalization  into  to  the  polyene. 

Table  IE.  Selected  bond  lengths  (r,  in  A,  according  to  labelling  scheme  below)  for  the 
neutral/charged  forms  of  1-3 


(CH3)2N+ 

o 

a 

CHs 

CHs 


c 

c 

c 

N 

N+ 


1(2+) 

2 

2(2.) 

3 

3(2+) 


1.45 

1.47 

1.45 

1.39 

1.37 


1.36 
1.34 

1.37 
1.36 
1.39 


1.44 
1.46 
1.41 

1.45 
1.41 


1.37 

1.36 

1.43 
1.36 

1.44 


1.43 

1.44 
1.35 
1.44 
1.35 


apirst  bond  adjacent  to  ring;  *>Second  bond  adjacent  to  ring. 


Molecule  2  is  quinoidal  in  the  neutral  state,  but  upon  reduction,  becomes  more  aromatic 
(Table  IE).  The  C-0  bond  (r(0-Csp2))  increases  from  1.24  A  to  1.26  A  on  reduction,  consistent 
with  a  lowering  of  bond  order  suggested  from  the  resonance  structures  in  Figure  1.  Furthermore, 
the  center  of  the  molecule  displays  the  geometry  expected  from  the  aromatic  form  of  2(2-)  with 
the  central  single  and  double  bond  lengths  of  that  form  at  1.44  A  and  1.35  A,  respectively. 
Conversion  to  the  aromatic  form  is  not  complete  since,  all  of  the  negative  charge  would  localize 
solely  upon  the  oxygen  atoms.  Nevertheless,  the  expected  decrease  in  a  and  y  due  to  the  increase 
in  aromaticity  should  partially  offset  increases  in  a  or  y  due  to  the  two  electron  reduction. 


Correspondingly,  much  smaller  increases  in  a  and  y  are  predicted  upon  reduction  of  2  than  are 


predicted  for  1  upon  oxidation.  The  increase  does,  however,  suggest  that,  for  this  molecule,  the 
charges  at  the  ends  of  the  molecule  and  the  decrease  in  bond  length  alternation  are  more  important 
than  the  increase  in  aromaticity.  When  3  is  oxidized,  its  geometry  distinctly  changes  from 
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quinoidal  to  aromatic  (Table  III)  and  both  oc  and  y  decrease  (Table  I).  Thus,  here,  oxidation 
(bipolaron  formation)  increases  both  aromaticity  and  to  a  smaU  extent  bond  length  alternation,  both 
of  which  appear  to  reduce  a  and  y. 

In  summary,  for  polyenes,  4  and  5,  bipolaron  (charge)  formation  leads  to  enhanced  a  and 
y.  In  the  ring  containing  compounds  there  are  two  effects:  the  charge  formation  which  may  be 
expected  to  increase  cc  and  y  and  a  change  from  aromatic  to  quinoidal  structures  (or  vice  versa) 
concomitant  with  bipolaron  formation.  The  electrons  in  the  molecules  with  aromatic  endgroups  are 
expected  to  be  less  delocalized,  resulting  in  a  detrimental  effect  on  y.  The  molecular  topology 
determines  whether  these  two  effects,  charge  formation  and  change  in  aromaticity  are  constructive 
or  competitive.  In  1,  upon  oxidation,  charges  are  formed  and  aromaticity  is  reduced; 
consequently,  the  largest  increase  in  y  is  predicted.  In  contrast  for  2  and  3  the  aromaticity 
increases  upon  charge  formation,  and  the  two  effects  cancel  to  a  certain  extent.  Thus,  in  2,  only  a 
small  increase  in  y  is  predicted  upon  reduction,  and  in  3  y  is  predicted  to  decrease  upon  oxidation. 
In  conclusion,  from  these  examples  we  suggest  that  it  is  not  just  the  oxidation  state  (and  charge)  of 
a  molecule  but  also  its  topology,  degree  of  aromaticity  and  bond  length  alternation  that  influence  a 

and  y. 


CYANINES 


The  fact  that  cyanines  exhibit  very  little  bond  length  alternation  can  be  rationalized  by 
considering  their  two  resonance  structures  (Figure  2).  Since  the  two  resonance  forms  are 
degenerate,  each  is  expected  to  make  an  equal  contribution  to  the  description  of  the  ground-state  of 
the  molecule,  leading  to  no  bond  length  alternation.  Short  chain  length  cyanines,  such  as  6  (n  = 
2  and  3)  have  large  negative  y  values.^®  Negative  y  is  found  in  other  related  molecules  with 

degenerate  resonance  structures  such  as  squarylium  dyes.^^ 


(CH3)2+n 


(CHalzN 


(CHalz 


Figure  2.  The  two  resonance  forms  contributing  to  the  structure  of  6,  shown  for  n  =  2. 

For  short  chain  length  cyanines,  (6),  our  computations  predict  that  y  is  negative  when  n  is 
small  and  increases  in  chain  length,  initially  increase  both  a  and  lyl.  The  correspondence  between 
increased  polarizabilities  and  increased  molecular  length  is  ubiquitous  in  NLO  molecular  structure- 
property  studies.  Interestingly,  dramatic  changes  in  a  and  y  occur  when  n  increases  from  6  to  7. 
The  sign  of  y  changes,  a  decreases,  and  the  bond  length  alternation  (Figure  3)  increases 
substantially.  From  the  geometry,  it  appears  that  charge  is  localizing  on  one  side  of  these  longer 
molecules,  indicating  that  the  two  resonance  structures  are  contributig  unequally  to  the  ground- 
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state  description  of  the  molecule.  The  extent  of  charge  delocalization  in  long  conjugated  systems 
has  been  considered  for  doped  conjugated  polymers^®  and  model  compounds. Using  other 
parameterizations  available  in  MOPAC  (MNDO,  PM3),  we  observed  increases  in  bond  length 
alternation  and  corresponding  changes  in  a  and  y  for  6  when  n  =  6  and  5,  respectively.  Thus,  this 
computation  makes  an  interesting,  specific  prediction  about  the  chain  length  at  which  cyanines  will 
start  to  display  charge  localization  and  resulting  increases  in  bond  length  alternation.  The 
calculations  suggest  that,  beyond  this  point,  y  becomes  increasingly  positive,  and  for  n  =  8  it  in 
fact  changes  sign.  Certainly,  the  point  at  which  these  changes  take  place  will  be  sensitive  to  the 
parameterization  of  the  semi-empirical  method  and  thus  these  results  should  be  viewed  with 
caution.  Thus,  it  will  be  of  interest  to  synthesize  long  cyanines  and  experimentally  pursue  this 
phenomenon  as  it  suggests  that  the  ultimate  volume  unit  values  for  a  and  lyl  is  limited.  Finally,  we 
note  that  at  comparable  lengths,  cyanines  exhibit  much  large  lyl  than  that  of  bipolarons.  Thus, 
although  6  (n  =  7)  is  one  atom  shorter  that  4  (2+),  lyl  for  the  former  compound  is  more  than 
four  times  that  of  the  latter. 


Figure  3.  Bond  length  alternation  (defined  as  the  difference  in  the  average  of  all  r(C-C)  distances 
minus  the  average  of  all  r(C=C)  distances  in  6  for  n  =  2-10  (Figure  1). 
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